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j~^ , Perturbation theory for gravity in dimensions greater than two requires 

Cn I higher derivatives in the free action. Higher derivatives seem to lead to ghosts, 

^vq I states with negative norm. We consider a fourth order scalar field theory and 

K^ ' show that the problem with ghosts arises because in the canonical treatment, 

QQ ■ (j) ^-iid 0(j) are regarded as two independent variables. Instead, we base quan- 

^^ , tum theory on a path integral, evaluated in Euclidean space and then Wick 

L». i rotated to Lorentzian space. The path integral requires that quantum states 

be specified by the values of cj) and (p^r- To calculate probabilities for obser- 
vations, one has to trace out over (p^r on the final surface. Hence one loses 
unitarity, but one can never produce a negative norm state or get a negative 
JL , probability. It is shown that transition probabilities tend toward those of the 

D i second order theory, as the coefficient of the fourth order term in the action 



tends to zero. Hence unitarity is restored at the low energies that now occur 
in the universe. 



I. INTRODUCTION 

In standard, second order theory the Lagrangian is a function of the fields and their 
first derivatives. The path integral is calculated by perturbation theory, with the part of 
the action that contains quadratic terms in the fields and their first derivatives regarded as 
the free field action, and the remaining terms as interactions. One then calculates Feynman 
diagrams, using the interactions as vertices, and the propagator defined by the free part 
of the action. This is equivalent to calculating the expectation value of the interactions in 
the Gaussian measure defined by the free action. One would therefore expect perturbation 
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theory to make sense, when and only when, the interaction action is bounded by the free 
action. 

This is born out by the examples we know. In two dimensions, the free action of a scalar 
field (t>. 



S= dx 



in0 + m'^d)'^ 



(1) 



is the first Sobolev normQ ll^lb.i of the field 0. In two dimensions, the first Sobolev norm 
bounds the pointwise value of 0, thus it also bounds the volume integral of any entire 
function of (p. This means that the free action bounds any interaction action, so perturbation 
theory should work. Indeed one finds that in two dimensions, any quantum field theory is 
renormalizable. 

In four dimensions on the other hand, the first Sobolev norm does not bound the point- 
wise value of (J), but only the volume integral of 0^. This means that the free action bounds 
the interactions only for theories with quartic interactions, like A0^, or Yang-Mills. Indeed, 
these are the quantum field theories that are renormalizable in four dimensions. Note that 
even Yang-Mills is not renormalizable in dimensions higher than four, because the interac- 
tions are not bounded by the free action. Similarly, Born-Infeld is not renormalizable in 
dimensions higher than two. 

When one does perturbation theory for gravity, one writes the metric as go + 6g, where 
qq is a background metric that is a solution of the field equations. The terms quadratic in Sg 
are again regarded as the free action, and the higher order terms are the interactions. The 
latter include terms like {'V6g)'^, multiplied by powers of 6g. The volume integral of such 
an interaction is not bounded by the free action and perturbation theory breaks down for 
gravity, which is not renormalizable . Even if all the higher loop divergences canceled by 
some miracle in a supergravity theory, one couldn't trust the results, because one is using 
perturbation theory beyond its limit of validity; 6g can be much larger than gQ locally for 
only a small free action. In other words, there are large metric fluctuations below the Planck 
scale. 

The situation is different however if one adds curvature squared terms to the Einstein- 
Hilbert action. The action is now quadratic in second derivatives of 6g, so one takes the free 
action to be the quadratic terms in Sg, and its first and second derivatives. This means that 
it is the second Sobolev norm 115(7112,2 of Sg, which bounds the pointwise value of Sg. Hence 
the free action bounds the interactions, and perturbation theory works. This is reflected in 
the fact that the R + B? theory is renormalizable , and in fact asymptotically free [Q . 



^For a function / G C'^{M), 1 < p < oo, and an integer A: > 0, the Sobolev norm is defined [Q] as 
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p,k 



'*^ 0<J<fc 



(2) 



where \D^ f\ is the pointwise norm of the jth covariant derivative and /ig is the Riemannian volume 
element. 



However, higher derivatives seem to lead to ghosts, states with negative norm, which have 
been thought to be a fatal flaw in any quantum fleld theory (see e.g. p). 

In the next section we review why higher derivatives appear to give rise to ghosts. The 
existence of ghosts would mean that the set of all states would not form a Hilbert space 
with a positive deflnite metric. There would not be a unitary S matrix, and there would 
apparently be states with negative probabilities. These seemed sufficient reasons to dismiss 
any quantum fleld theory, such as Einstein gravity, that had higher derivative quantum 
corrections and ghosts. However, we shall show that one can still make sense of higher 
derivative theories, as a set of rules for calculating probabilities for observations. But one 
can not prepare a system in a state with a negative norm, nor can one resolve a state into 
its positive and negative norm components. So there are no negative probabilities, and no 
non unitary S matrix. 

Although gravity is the physically interesting case, in this paper we consider a fourth 
order scalar fleld theory, which has the same ghostly behaviour, but doesn't have the com- 
plications of indices or gauge invariance. We show explicitely that the higher derivative 
theory tends toward the second order theory, as the coefficient of the fourth order term in 
the action tends to zero. Hence the departures from unitarity for higher derivative gravity 
are very small at the low energies that now occur in the universe. 

II. HIGHER DERIVATIVE GHOSTS 

We consider a scalar fleld with a fourth-order Lagrangian in Lorentzian signature, 

L = -U (n - m^) (n - m^) - A0^ (3) 



where m2 > rrii. Deflning 



[n-ml)(t> (n-ml) 



the Lagrangian can be rewritten as 

L = \lj, (n - m?) ^1 - \^, (n - mt) ^2 - ^^/^^,^, (V^i - ^2Y (5) 

The action of ■02 has the wrong sign. Classically it means that the energy of the ■02 
fleld is negative, while that of "01 is positive. If there were no interaction term, this negative 
energy wouldn't matter because each of the flelds, V^i and 02? would live in its own world and 
the two worlds would not communicate with each other. However, if there is an interaction 
term, like 0"^, it will couple i\}\ and 02 together. Energy can then flow from one to the other, 
and one can have runaway solutions, with the positive energy of V'l and the negative energy 
of 0^2 both increasing exponentially. 

In quantum theory, on the other hand, one is in trouble even in the absence of interac- 
tions, as can be seen by looking at the free field propagator for 0. In momentum space, this 
is the inverse of a fourth order expression in p, which can be expanded as 



{ml — ml) \{p'^ + nil) (p^ + ^^i)/ ' 

This is just the difference of the propagators for ipi and 1^2- The important point is that 
the propagator for 1^2 appears with a negative sign. This would mean that states with an 
odd number of 1^2 particles, would have a negative norm. In other words, ip2 particles are 
ghosts. There wouldn't be a positive definite Hilbert space metric, nor a unitary S matrix. 

If there weren't any interactions, the situation wouldn't be too serious. The state space 
would be the direct sum of two Hilbert spaces, one with positive definite metric and the 
other negative. There wouldn't be any physically realized operators that connected the 
two Hilbert spaces, so ghost number would be conserved by a superselection rule. A 0^ 
interaction however, would allow tp2 particles to be created or destroyed. As in the classical 
theory, there will be instabilities, with runaway production of V'l and tp2 particles. These 
instabilities show up in the fact that interactions tend to shift the ghost poles in the two 
point function for (p into the complex p-plane, where they represent exponentially growing 
and decaying modes |§J31- 

It seems to add up to a pretty damning indictment of higher derivative theories in general, 
and quantum gravity and quantum supergravity in particular. However, the problem with 
ghosts arises because in the canonical treatment, and D^ are regarded as two independent 
variables, although they are both determined by 0. We shall show that, by basing quantum 
theory on a path integral over the field, evaluated in Euclidean space and then Wick rotated 
to Lorentzian space, one can obtain a sensible set of rules for calculating probabilities for 
observations in higher derivative theories. 

III. EUCLIDEAN PATH INTEGRAL 

According to the canonical approach, one would perform the path integral over all ipi 
and tp2- The path integral over tpi will converge, but the path integral over ?/'2 is divergent, 
because the free action for 1JJ2 is negative definite. However, one shouldn't do the path 
integrals over ipi and -02 separately because they are not independent fields, they are both 
determined by 0. The fourth order free action for is positive definite, thus the path integral 
over all in Euclidean space should converge, and should define a well determined Euclidean 
quantum field theory. 

One way to compute the path integral for a fourth order theory, is to expand in 
eigenfunctions of the differential operator O in the action. One then integrates over the 
coefficients in the harmonic expansion, which gives (detO)~^/^. Another way is to use time 
slicing, by dividing the period into a number of short time steps e and approximating the 
derivatives by 

, (0n+l - 0n) , (0n+2 - 20n+i + 0n) 

e e^ 

One then integrates over the values of on each time slice. In a second order theory, where 
the action depends on and 0^^ but not on (p^rr, the path integral will depend on the values 
of on the initial and final surfaces. However, in a fourth order theory, the use of three 



neighbor differences means that one has to specify (p^r on the initial and final surfaces as 
well. 

One can also see what needs to be specified on the initial and final surfaces as follows. 
In classical second order theory, a state can be defined by its Cauchy data on a spacehke 
surface, i.e. the values of (p and (p^r on the surface. In a canonical 3+1 treatment, these 
are regarded as the position of the field and its conjugate momentum. In quantum theory, 
position and momentum don't commute, so instead one describes a state by a wave function 
in either position space or momentum space. In ordinary quantum mechanics, the position 
and momentum representations are regarded as equivalent: one is just the Fourier transform 
of the other. However, with path integrals, one has to use wave functions in the position 
representation. This can be seen as follows. Imagine using the path integral to go from a 
state at ti to a state at T2, and then to a state at t^. In the position representation, the 
amplitude to go from a field cpi on ri, to 02 at T2, is given by a path integral over all fields 
with the given boundary values. Similarly, the amplitude to go from 02 at T2, to 03 at t^, 
is given by another path integral. These amplitudes obey a composition law, 

G(03, 0l) = J ^02^(03, 02)^(02, 0l) (8) 

The composition law holds, only because one can join a field from 0i to 02 to a field 
from 02 to 03, to obtain a field from 0i to 03. Although in general 0^,- will be discontinuous 
at t2, the field will still have a well defined action, 

5(03,00= 5(03, 02) +5(02, 03) (9) 

On the other hand, if one would use the momentum representation and wave functions in 
terms of 0^,-, the composition law would no longer hold, because the discontinuity of at T2 
would make the action infinite. Thus in second order theories, one should use wave functions 
in terms of rather than 0^^. 

In a fourth order theory, a classical state is determined by the values of and its first 
three time derivatives on a spacelike surface. In a canonical treatment, and 0t-t are 
usually taken to be independent coordinates. For the scalar field theory (|^) we then have 
the conjugate momenta 

This suggests that in quantum theory, one should describe a state by a wave functional 
^(0,0,rT) on a surface. Indeed, this is closely related to using the fields ipi and -02 that 
we introduced earlier. These were linear combinations of and ^0, thus taking the wave 
function to depend on ipi and -02, is equivalent to it depending on and (p^rr- However, 
if one does the path integral between fixed values of and (p^rr, one gets in trouble with 
the composition law, because the values of (p^r on the intermediate surface at T2 are not 
constrained. Hence (p^r will be in general discontinuous at T2, which implies that (p^rr will 
have a delta-function when one joins the fields above and below T2. In a second order action 
,-r appears linearly, thus the delta-function can be integrated by parts and the action of the 
combined field is finite. But in a fourth order action {(p^rr)^ appears, rendering the action 
of the combined field infinite if rr is a delta-function. 



Therefore, the path integral requires that quantum states be specified by and (p^r in 
order to get the composition law for amplitudes in a fourth order theory. In the next section 
we show how one can obtain transition probabilities for observations from the Euclidean 
path integral over (p. 

IV. HIGHER DERIVATIVE HARMONIC OSCILLATOR 

A. Ground State Wave Function 

To illustrate how probabilities can be calculated, we consider a higher derivative harmonic 
oscillator, for which in Euclidean signature we take the action 



S 



dr 



a" ,2 i 



-m 
2 



2 J,2 



(11) 



For a^ > 0, this is very similar to our scalar field model, since in the latter we can take 
Fourier components so that spatial derivatives behave like masses. The general solution to 
the equation of motion is given by 



(r) = A sinh Air + B cosh Air + C sinh A2r + D cosh A2r, 



(12) 



where Ai and A2 are given by (|A2|) . For small a, Ai ~ -m and A2 ~ 1/a. 

The fourth order action for (p is positive definite, thus it gives a well defined Euclidean 
quantum field theory. In this theory, one can calculate the amplitude to go from a state 
(0i)0i,r) at time ri, to a state (025 02,r) at time r2. In particular, one can calculate the 
ground state wave function, the amplitude to go from zero field in the infinite Euclidean 
past, up to the given values (0o, 0o,r) at r = 0. This yields (see Appendix A) 



^ 



o{(Po,(Po,t) 



where 



iV'exp 



F' 



-F' 



^0,T 



+ 



m 



a 



+ 



2m^ 



(A2 - Ai) 



m/a 

— ^0O0O,r 



:i - Am^a^] 



2«2(Ai + A2)(A2-Ai 



(13) 



(14) 



and N'{a,m) is a normalization factor. 

Similarly, one can calculate the Euclidean conjugate ground state wave function \E'q, the 
amplitude to go from the given values at r = 0, to zero field in the infinite Euclidean future. 
This conjugate wave function is equal to the original ground state wave function, with the 
opposite sign of 0o,t- The probabihty that a quantum fiuctuation in the ground state gives 
the specified values 0o and 0o r on the surface r = 0, is then given by 



P( 



^0, (PO,r) 



^0^* = N" exp 



-2F' 



^l^'^^l 



(15) 



The probability dies off at large values of (p and (p^r and is normalizable, thus the prob- 
ability distribution in the Euclidean theory is well-defined. However if one Wick rotates 
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to Minkowski space, 0^^ picks up a minus sign. The probability distribution becomes un- 
bounded for large Lorentzian .^ and can no longer be normalized. This is another reflection 
of the same problem as the ghosts. You can't fully determine a state on a spacelike surface, 
because that would involve specifying and Lorentzian 1, which doesn't have a physically 
reasonable probability distribution. 

Although one can not deflne a probability distribution for and Lorentzian 0^t on a 
spacelike surface, one can calculate a probability distribution for alone, by integrating out 
over Euclidean ,-. This integral converges because the probability distribution is damped 
at large values of Euclidean 0,-. This is just what one would calculate in a second order 
theory. So the moral is, a fourth order theory can make sense in Lorentzian space, if you 
treat it like a second order theory. The normalized probability distribution that a ground 
state fluctuations gives the specifled value 0o on a spacelike surface is then given by. 



Pi 



'2F'm^ 



-na 



1/2 



exp 



2mF' 



a 



(16) 



As the coefficient a of the fourth order term in the action tends to zero, this becomes 



m 



1/2 



ma. 



7J (1+ 2 



exp[— m(l + ma)0o]. 



(17) 



which tends toward the result for the second order theory. 



B. Transition Probabilities 



In this section we compute the Euclidean transition probability, to go from a specified 
value 01 at time ri, to 02 at time r2, for the higher derivative harmonic oscillator. 

In a second order theory, a state can be described by a wave function that depends on the 
values of on a spacelike surface. Thus a transition amplitude is given by a path integral 
from an initial state 0i on ri, to a final state 02 on T2- To calculate the probability to go 
from the initial state to the final, one multiplies the amplitude by its Euclidean conjugate. 
This can be represented as the path integral from a third surface, at ts, back to T2. Because 
the path integral in a second order theory depends only on on the boundary, what happens 
above r^ and below Ti doesn't matter. Furthermore, the path integrals above and below T2 
can be calculated independently, which implies the probability to go from initial to final, 
can be factorized into the product of an S matrix and its adjoint. The S matrix is unitary, 
because probabihty is conserved. 

Now let us calculate the probability to go from an initial to a final state in the fourth 
order theory ([TT|). The path integral requires quantum states to be specified by and 0^^. 
The transition amplitude to go from a state (0i, 01,,-) at time ri = — T, to a state (02, 02,t) 
at time r2 = 0, reads 

((02, 02,.; 0)1(01, 01,.; -T)) = /^'^''^'"^rf[0(r)]exp[-5(0)] (18) 

This is evaluated in Appendix A, by writing = 0d + 0', where 0c/ obeys the equation of 
motion with the given boundary conditions on both surfaces. 



The result is 



((</>2,<^2,r;0)|(</.i,</.i,.;-T)) 



a{l + aN)H' 
2^ 



1/2 



exp 



-E{4>1 + <pl) - F{4>1 + <pU 



- ^01,^02,^ + H(j)i(j)2 - K{(f)2,T4>2 - 4>1,t4>i) - L{4>2,t4>1 " 4>1,t4>2)] (19) 

The coefficient functions in the exponent are given by (|A6|) , and A^ is a normahzation factor. 

Again, one can construct a three layer 'sandwich' to calculate the probability to go from 
the initial state to the final. However, in contrast with the second order theory the path 
integral now depends on both (p and (j),r on the boundaries. This has two important im- 
plications for the calculation of the transition probability. Firstly, as we just showed, one 
can't observe Lorentzian (p^j. because it has an unbounded Lorentzian probability distribu- 
tion. Therefore one should take (p^r to be continuous on the surfaces and integrate over all 
values, fixing only the values of (p on the surfaces. Because the path integrals above and 
below T2 = both depend on 4>2,t, the probability P{4>2, 4>i) to observe the initial and final 
specified values of (p does not factorize into an S matrix and its adjoint. Instead, there is loss 
of quantum coherence, because one can not observe all the information that characterizes 
the final state. 

After multiplying by the Euclidean conjugate amplitude and integrating out over 02,t we 
obtain 



a{l + aN)H 



27r2 



n \V2 
2f) ^"P 



-2E( 



^\ + 



2F(Pl^ + 2H(Pv 



')2-\ <; 

2F^ 



ir 



(20) 



Another consequence of the dependence of the path integral on (p^j. is that what goes on 
outside the sandwich, now affects the result. The most natural choice, would be the vacuum 
state above t^ = T and below ti = —T. In other words, one takes the path integral to be over 
all fields that have the given values on the three surfaces, and that go to zero in the infinite 
Euclidean future and past. This means that to obtain the transition probability we also 
ought to multiply by the appropriately normalized ground state wave function \l/o(0i, 0i,t) 
and its Euclidean conjugate. The probability P{(p2, <Pi) is then given by 



P( 



^2,(Pl 



d[0i,,]^O^Oy C?[02,r]((01, 01,r)|(02, 02,r))((02, 02,r)|(01, 01,r)) 



a'il 



27r2(4F(F' + F) - G2) 



exp 



-2Ei(Pl + (Pl) 



TflF 

2 (Pi + 2i70i02 



a 



Here F'{a, m) is the coefficient in the exponent of the ground wave function 
a normalization factor. In the limit a —>■ 0, this reduces to 



771 



^2,(Pl, 



27r sinh mT 



exp 



mcosh.mT{(p1 + 



sinh mT 



2m(pi(p2 ,2 
mcp]^ 



(21) 
and A^ is 

(22) 



Hence the probability given by the sandwich tends toward that of the second order 
theory, as the coefficient of the fourth order term in the action tends to zero. This is 
important, because it means that fourth order corrections to graviton scattering can be 
neglected completely at the low energies that now occur in the universe. On the other hand, 
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in the very early universe, when fourth order terms are important, we expect the Euchdean 
metric to be some instanton, hke a four sphere. In such a situation, one can not define 
scattering or ask about unitarity. The only quantities we have any chance of observing are 
the n-point functions of the metric perturbations, which determine the n-point functions of 
fluctuations in the microwave background. With Reall we have shown that Starobinsky's 
model of inflation [§], in which inflation is driven by the trace anomaly of a large number of 
conformally coupled matter fields, can give a sensible spectrum of microwave fluctuations, 
despite the fact it has fourth order terms and ghosts 0. Moreover, the fourth order terms 
can play an important role in reducing the fluctuations to the level we observe. 

Finally, in order to obtain the Minkowski space probability, one analytically continues T2 
to future infinity in Minkowski space, and ri and ts to past infinity, keeping their Euclidean 
time values fixed. This gives the Minkowski space probability, to go from an initial value 0i 
to a final value (ho. 



V. RUNAWAYS AND CAUSALITY 

The discussion in Section II suggests that even the slightest amount of a fourth order 
term will lead to runaway production of positive and negative energies, or of real and ghost 
particles. The classical theory is certainly unstable, if one prescribes the initial value of and 
its first three time derivatives. However, in quantum theory every sensible question can be 
posed in terms of vacuum to vacuum amplitudes. These can be defined by Wick rotating to 
Euclidean space and doing a path integral over all fields that die off in the Euclidean future 
and past. Thus the Euclidean formulation of a quantum field theory implicitly imposes 
the final boundary condition that the fields remain bounded. This removes the instabilities 
and runaways, like a final boundary condition removes the runaway solution of the classical 
radiation reaction force. The price one pays for removing runaways with a final boundary 
condition, is a slight violation of causality. For instance, with the classical radiation reaction 
force, a particle would start to accelerate before a wave hit it. This can be seen by considering 



a single electron which is acted upon by a delta-function pulse |]10[. The equation of motion 
for the x-component reduces to 

xtt = Ax,iti + 5(t), (23) 

with A = S^- This has the solution 

,(t) = l^eM-,^t]-^,. (24) 

The integrand has two singularities, at cj = and uj = iX^^. The final boundary condition 
that X t should tend to a finite limit, implies one must choose an integration contour that 
stays close to the real axis, going below the second singularity. This yields 

x(t) = Aexp[t/A], t<0 

= t + X, t > (25) 

which is without runaways, but acausal. 



However, this pre-acceleration is appreciable only for a period of time comparable with 
the time for light to travel the classical radius of the electron, and thus practically unob- 
servable. 

Similarly, if we would add an interaction term to the higher derivative scalar field theory 
(0), the imposition of a final boundary condition to eliminate the runaway solutions, would 
lead to acausal behaviour on the scale of m^^, where m2 is the mass of the ghost particle. 
However, in the context of quantum gravity, one could again never detect a violation of 
causality, because the presence of a mass introduces a logarithmic time delay At ~ —m log 6, 
where h is the impact parameter. Thus there is no standard arrival time, one can always 
arrive before any given light ray by taking a path which stays a sufficiently large distance 
from the mass. 

VI. CONCLUDING REMARKS 

We conclude that quantum gravity with fourth order corrections can make sense, despite 
apparently having negative energy solutions and ghosts. In doing this, we seem to go against 
the convictions of the last 25 years, that unitarity and causality are essential requirements of 
any viable theory of quantum gravity. Perturbative string theory has unitarity and causal- 
ity, so it has been claimed as the only viable quantum theory of gravity. But the string 
perturbation expansion does not converge, and string theory has to be augmented by non 
perturbative objects, like D-branes. One can have a world-sheet theory of strings without 
higher derivatives, only because two dimensional metrics are conformally flat, meaning per- 
turbations don't change the light cone. Still, we live either on a 3-brane, or in the bulk of a 
higher dimensional compactifled space. The world-sheet theory of D-branes with p greater 
than one has similar non-renormalizability problems to Einstein gravity and supergravity. 
Thus string theory effectively has ghosts, though this awkward fact is quietly glided over. 

To summarize, we showed that perturbation theory for gravity in dimensions greater 
than two required higher derivatives in the free action. Higher derivatives seemed to lead 
to ghosts, states with negative norm. To analyze what was happening, we considered a 
fourth order scalar fleld theory. We showed that the problem with ghosts arises because in 
the canonical approach, and 00 are regarded as two independent coordinates. Instead, 
we based quantum theory on a path integral over 0, evaluated in Euclidean space and then 
Wick rotated to Lorentzian space. We showed the path integral required that quantum 
states be specifled by the values of and 0^^ on a spacelike surface, rather than and 
0_^T as is usually done in a canonical treatment. The wave function in terms of and 0^,- is 
bounded in Euclidean space, but grows exponentially with Minkowski space 0,t-- This means 
one can not observe 0^^ but only 0. To calculate probabilities for observations one therefore 
has to trace out over 0_^ on the flnal surface, and lose information about the quantum state. 
One might worry that integrating out 0^^ would break Lorentz invariance. However, ,- is 
conjugate to (p^rr so tracing over 0^^ is equivalent to not observing 00. Since, according to 
eq.(^), '01 and V'2 are linear combinations of and ^0, this means that one only considers 
Feynman diagrams whose external legs are ipi — ip2- You don't observe the other linear 
combination, m\il)i — m\%l)2. 

Because one is throwing away information, one gets a density matrix for the flnal state, 
and loses unitarity. However, one can never produce a negative norm state or get a negative 
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probability. We illustrated with the example of a higher derivative harmonic oscillator that 
probabilities for observations tend toward those of the second order theory, as the coefficient 
of the fourth order term in the action tends to zero. This means that the departures from 
unitarity for higher derivative gravity will be very small at the low energies that now occur 
in the universe. On the other hand, the higher derivative terms will be important in the 
early universe, but there unitarity can not be defined. 

Acknowledgements 

It is a pleasure to thank David Gross, Jim Hartle and Edward Witten for helpful dis- 
cussions. We would also like to thank the ITP at Santa Barbara and the Department of 
Physics at Caltech, where some of this work was done, for their hospitality. TH is Aspirant 
FWO, Belgium. 

APPENDIX A: TRANSITION AMPLITUDE 

We compute the Euclidean transition amplitude, to go from an initial state (0i, 0i_^) on 
a spacelike surface at r = — T, to a final state (02,02,t) at r = 0, for the higher derivative 
harmonic oscillator (pTI). The general solution to the equation of motion is given by 

0(r) = A sinh Air + B cosh Air + C sinh A2r + D cosh A2r, (Al) 

where 

Ai = ^i=V(l-Vl-4mV) , X2 = ^=\/{l + Vl-4m^a^) (A2) 

v2a"' v2a^ 

The transition amplitude is given by a path integral, 

/i'f>2,<t>2,T) 
d[<p{T)]exp[-Si<P)] (A3) 

This can be evaluated by separating out the 'classical' part of (p. If we write (p = (f>ci + </>', 
where (pci obeys the equation of motion with the required boundary conditions on both 
surfaces r = and t = T, then the amplitude becomes 

((02,02,.;O)|(0i,0i,.;-r)) = exp[-^,K0i,0i,r,02,02,.)] r'^ rf[0'(r)]exp[-^(0')] (A4) 

J(0-T) 



The classical action is 

rT 



/ 2 / 2 2 y 2 



Sci = / dr 
Jo 

(.2 I j2\ I 77/i2 , ,2 



E{4>i + 02) + ^(0lr + 02,r) + ^01,^02,^ - i^0102 

+ -^(02,r02 - 01,T0l) + -^(02,T01 - 01,r02) (A5) 



where 
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E 



(cosh X2T — cosh AiT) (Ai sinh AiT + A2 sinh A2T) 



— m(l — Arri^Q^) 

2a^{\l - Af)P2 

1 1 

+ sinh AiT sinh A2T(Ai (2A^ + — ) cosh A2T sinh \{T - A2(2A^ + — ) sinh A2r cosh AiT) 



2m 



a 



a" 



a" 



(cosh A2r — cosh AiT) (A2 sinh AiT + Ai sinh A2T) 



(1 - Am^a^) 
2a\\l - Af)P2 

+ sinh AiT sinh A2T(A2(2A^ H — 5-) cosh A2T sinh AiT — Ai(2A2 H — 5-) sinh A2T cosh AiT) 



2m 



a 



a" 



a" 



G 



1 — Am^a 



2^,2n 



a 



2(Ai - A?)P2 V 



I — (cosh A2r cosh AiT — 1) 



a 



sinh AiT sinh A2T ) ( Ai sinh A2T — A2 sinh AiT) 



a^ 



H 



-"^^^-'"^'"'U^(coshA2TcoshA,T-l) 



a\\l - Af)p2 V a 



sinh AiT sinh A2T j (Ai sinh AiT — A2 sinh A2T) 



a^ 



K 



P2 



m / . 2 1 



a 



a" 



4m H — - I sinh AiT sinh A2T(1 — cosh Ai cosh A2T) 



2m2 



a^ 



(2 - 3(cosh2 AiT + cosh^ A2r)) 



-m(l-4mV) ^2m 
^ = a3(Ai - Af )P^ (^— (coshA2TcoshAiT - 1) 



sinh AiTsinh A2r j (cosh A2r — cosh AiT) 



with 



(A6) 

P ={\\ + \l) sinh AiT sinh A2T + 2AiA2(l - cosh AiT cosh A2T) (A7) 

The pre-exponential factor in ( [A4D can be derived from the classical action alone |TI| , it is 
basically the Jacobian of the change of variables (tti, (pi) — >■ ((/)2, (pi). Because the Lagrangian 
is quadratic, the prefactor is independent of the values specifying the initial and final states, 
and the transition amplitude ([A4| ) is exact. It is given by 



((02, 02,r; 0)1(01, 01,.; -T)) 



27r2 



exp 



-E{(Pi + (Pi) - P(0f , + 01J 



- G01,r02,T + -f^0102 - -f^(02,r02 " 01,t0i) " L{4>2,t4>1 " 01,r02)] (A8) 

The normalization factor A^ is independent of a to first order. It is determined by 
taking T -^ +cx) in (^) and requiring that the amplitude tends toward the product of two 
normalized ground state wave functions \I'o(0i, 0i,t) and \E'o(02, 02,t)- 

For small a, Ai ~ m and A2 ~ l/«, hence the transition amplitude becomes 



((02, 02,r; 0)1(01, 01,.; -T)) 



m,a 



277^ sinh TJiT 



1/2 



exp 



mcoshmT(0^ + 



a 



2 sinh m^T 



Xr ^ V^2,tJ 



ma COshmT(02,r02 - 01,r0l) m{a(p2,T + 02)(a01,r - 01 



sinh m,T 



sinh m^T 



(A9) 
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